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Abstract 

We first state a special type of Ito formula involving stochastic integrals of both standard and frac- 
tional Brownian motions. Then we use Doss-Sussman transformation to establish the link between back- 
ward doubly stochastic differential equations, driven by both standard and fractional Brownian motions, 
and backward stochastic differential equations, driven only by standard Brownian motions. Following 
the same technique, we further study associated nonlinear stochastic partial differential equations driven 
by fractional Brownian motions and partial differential equations with stochastic coefficients. 
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1 Introduction 

The theory of backward stochastic differential equations (BSDEs) and that of fractional Brownian motion 
(fBm) had developed simultaneously in their own separated directions for many years until Bender [2] gave 
an explicit solution to a linear BSDE driven by fBm in 2005. In 2009 Hu and Peng [10] stated a more 
general theory on fractional BSDEs by using the so-called quasi-expectation, but their case is still limited. 
The non-semimartingale property of the fBm (except the case of Hurst parameter H = 1/2, where it becomes 
a Brownian motion) makes it thorny to handle. Being not a semimartingale means there is no martingale 
representation theory for the fBm, which is crucial in the general BSDE theory (see the pioneering work of 
Pardoux and Peng [15]). In Jing and Leon [11], we tried to combine the fBm and BSDEs in another way: 
We transformed a semilinear backward doubly stochastic differential equation (BDSDE) driven by both a 
standard and a fractional Brownian motions with H S (0, 1/2) into a BSDE without integral with respect to 
the fBm, which turns out easier to deal with. The integral w.r.t. the fBm in the BDSDE was interpreted in 
the sense of the extended divergence operator and the integral was supposed to be linear w.r.t. the solution 
process. This allowed to apply the efficient tool of nonanticipating Girsanov transformation, developed in 
Buckdahn [4]. However, this method is restricted to semilinear BDSDEs. 

In this paper, we deal with BDSDEs for which the integrand of the integral with respect to the fBm is 
not necessarily linear with the solution process, and the Hurst coefficient H is supposed to belong t o the 
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interval(l/2, 1). Unlike the more irregular case H < 1/2, the stochastic integrals with respect to an fBm 
with H > 1/2 can be defined in different ways. So they can be defined with the help of the divergence 
operator in the frame of the Malliavin calculus, see Decreusefont and Ustiinel [8] and Alos et al. [1] (Notice 
that the Wick-Ito integral defined in Duncan et al. [9] coincides with the first one). They can also be defined 
pathwise as generalized Riemann-Stieltjes integral (see Zahle [24] and [25]) or with the help of the rough 
path theory (see Coutin and Qian [7] ) . For a complete list of references we refer to the two books by Biagini 
et al. [3] and Mishura [13]. 

Our approach to BDSDEs with an fBM is inspired by the work of Buckdahn and Ma [5] . In their study 
of stochastic PDEs driven by a Brownian motion B the authors of [5] used BDSDEs driven by B as well 
as an independent Brownian motion; the integral with respect to B is interpreted in Stratonovich sense. 
This allowed the application of the Doss-Sussman transformation in order to transform the BDSDE into a 
BSDE without integral with respect to B. On the other hand, the pathwise integral with respect to the fBm 
plays a role which is comparable with that of the Stratonovich integral in the classical theory. Nualart and 
Ra§canu [14] used the pathwise integral to solve (forward) stochastic differential equations driven by an fBm. 
For some technical reasons (such as the lack of Holder continuity, see Remark 4.12), we shall make use of 
the Russo-Vallois integral developed by Russo and Vallois in a series of papers ([21], [21], [21], etc.). Under 
standard assumptions which allow to apply the Doss-Sussman transformation, we associate the BDSDE 
driven by both a standard Brownian motion W and an fBm B, 



U t s > x = $(X t > x )+ [ S f(r,X t r ' x ,U t r ' x ,V r t ' x )dr+ f g{U^)AB r - f V* I dW r , se^t], 1 (1.1 
Jo Jo Jo 



) 

Jo 

with the BSDE driven only by the Brownian motion W , 

= <D(X*' X ) + f ' fax?* X' x > I" Z^idWr, se[0,t}. (1.2) 

Jo Jo 

Here / will be specified in Section 4; it is a driver with quadratic growth in z. We point out that the 
classical BDSDEs were first studied by Pardoux and Peng [16] and our BSDE (1.2) is a quadratic growth 
BSDE, which was studied first by Kobylanski [12]. In the works of Pardoux and Peng [16] and Buckdahn 
and Ma [5], i.e., when the Hurst parameter H — 1/2, one can solve the BDSDE directly and get the square 
intcgrability of the solution process. However, in the fractional case (H ^ 1/2), to our best knowledge, there 
does not exist a direct way to solve the BDSDE (1.1), and as it turns out in Theorem 4.7, we can only get 
that the conditional expectation of J* \ Z t r > x \ 2 dr is bounded by an a.s. finite process. This is also the reason 
that instead of using the space of square integrable processes, we use the space of a.s. conditionally square 
intcgrable processes (see the definition of the space TL'^(M. d ) in Section 2). 

A celebrated contribution of the BSDE theory consists in giving a form of probabilistic interpretation, 
nonlinear Feynman-Kac formula, to the solutions of PDEs (see, for instance, Peng [18], Pardoux and Peng 
[17]). As indicated in Kobylanski [12], the quadratic growth BSDE (1.2) is connected to the semilinear 
parabolic PDE 

%(t,x) = J?u(t,x)- f(t,x,u(t,x),a(x) T £u(t lX )), (t,x) e (0,T) xl»; 
u(0,x) = <P(x), xeW 1 , ( ' 

where ££ is the infinitesimal operator of a Markov process. Hence, it is natural for us to consider the form 
of equation (1.3) after Doss-Sussman transformation and we prove that it becomes the following semilinear 
SPDE 

du{t,x)= [&u(t,x) - f{t,x,u(t,x),a{x) T £u{t,x))] dt + g(u(t,x))dB t , (t, x) G (0, T) x W 1 ; 
u(0,x) = $(x), xeR n . { ' 

We emphasize that this paper can not be considered as a generalization of [11]. The reason is that, firstly, 
the Hurst parameters are distinct; secondly, the stochastic integrals with respect to the fBm are of different 
types. 



1 f ■ I dW s indicates that the integral is considered as the Ito backward one. 
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We organize the paper as follows. In section 2, we recall some basic facts about the fractional Brownian 
motion, we give the general framework of our work, and we recall the definition of the backward Russo-Vallois 
integral as well as some of its properties. In section 3 we prove a type of ltd formula involving integrals with 
respect to both standard and fractional Brownian motions, which will play an important role in the following 
sections. We perform a Doss-Sussman transformation in Section 4 to transform a nonlinear BDSDE (1.1) 
into a BSDE (1.2) and show the relationship between their solutions. In particular, we show that BSDE 
(1.2) has a unique solution (Y t,x , Z t,x ), and the couple of processes (£7*> x , V t:X ) associated with (Y l ' x , Z l ' x ) 
by the inverse Doss-Sussman transformation is the unique solution of BDSDE (1.1). Finally, the stochastic 
PDE associated with BDSDE (1.1) is briefly discussed in Section 5. 



2 Preliminaries 

2.1 Fractional Brownian Motion and General Setting 

In this subsection we recall some basic results on the fBm and the related setting. For a more complete 
overview of the theory of fBm, we refer the reader to Biagini et al. [3] and Mishura [13]. 

Let (f2',J 7 ',P / ) be a classical Wiener space with time horizon T > 0, i.e., SI' = Co([0,T];R) denotes 
the set of real-valued continuous functions starting from zero at time zero, endowed with the topology of 
the uniform convergence, B(Sl') is the Borel cr-algebra on SI' and P' is the unique probability measure on 
(Sl,B(Sl')) with respect to which the coordinate process W t °(o/) = co'(t), t G [0, T], w' G SI' is a standard 
Brownian motion. By T' we denote the completion of B(SV) by all P'-null sets in SI' . Given H G (1/2, 1), 
we define 

B t = f K H (t,s)dW°, te [0,T], 
Jo 

where Kh is the kernel of the fBm with parameter H G (1/2, 1): 



K H (t, s) = Cms 1 / 2 -" f u H - l / 2 (u - S ) H ~ 3 / 2 du 

J s 



with Cr = y (2H-i)/3(2 H 2H H-i/2) • ^ * s wcn known that such defined process B is a one-dimensional fBm, 
i.e., it is a Gaussian process with zero mean and covariance function 

R H (t, s) := E [B t B s ] - \ (t 2H + s 2H -\t- s\ 2H ) , s,t e [0, T]. 

We let {Wt : < t < T} be the coordinate process on the classical Wiener space (f2", J 7 " ,P") with 
Q" = Co([0, T]; M d ), which is a d-dimensional Brownian motion with respect to the Wiener measure P". We 
put (SI, J" , P) = (SI', I', P') ® (SI", T", P") and let T = J" V A/", where TV is the class of the P-null sets. We 
denote again by B and W the canonical extensions of the fBm B and of the Brownian motion W from SI' 
and SI", respectively, to SI. 

We let F™ T] = a{W T -W s ,t<s<T} V7V, T? = a{B s , < s < t} VAf, and Q t = T^ T] V Tf , t e [0, T]. 
Let us point out that Fu^t] 18 decreasing and Tf is increasing in t, but Gt is neither decreasing nor increasing. 
We denote the family of cr-fields {Gt}o<t<T by G. Moreover, we also introduce the backward nitrations 
H = {Ht = V^} teM and ¥ w = {^ T] }te%T]- 

Finally, we denote by C(M, [0,r];M m ) the space of the R m -valued continuous processes {ift,t G [0,T]} 
such that ip t is "H t -measurable, t G [0,T], and M 2 (¥ w , [0, T}; R m ) the space of the Revalued square- 
integrable processes {ipt,t G [0,T]} such that ip t is J 7 ^] -measurable, t G [0,T]. Let (R) be the set of 
H-progressively measurable processes which are almost surely bounded by some real- valued J^-measurable 
random variable, and let H 2 ^(R d ) denote the set of all Revalued H-progressively measurable processes 



7 = {it ■ t G [0, T}} such that E J Q T | Tt | 2 dt|J^ 



< +oo, P-a.s. 
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2.2 Russo-Vallois Integral 

In a scries of papers ([20], [21], [22], etc.), Russo and Vallois denned new types of stochastic integrals, namely 
forward, backward and symmetric integrals, which are extensions of the classical Riemann-Sticltjes integral, 
and in fact these three integrals coincide, when the integrator is a fBm with Hurst parameter H € (1/2, 1). 
Here we will mainly use the backward Russo-Vallois integral in this paper. It turns out to be a convenient 
definition for stochastic integral with respect to our fBm B. 

Let us recall some results by Russo and Vallois which we will use later. In what follows, we make the 
convention that all continuous processes {X t , t € [0, T]} are extended to the whole line by putting X t = X , 
for t < 0, and X t = X T , for t > T. 

Definition 2.1. Let X and Y be two continuous processes. For e > 0, we set 

1 /"* 

I{e,t,X,AY)^ - / X(s)(Y(s)-Y(s-s))ds, 
£ Jo 

C £ (X,Y)(t)±- f (X(s)-X(s-e))(Y(s)-Y(s-s))ds, t e [0,T]. 
£ Jo 

Then the backward Russo-Vallois integral is defined as the uniform limit in probability as e + , if the limit 
exists. The generalized bracket [X, Y] is the uniform limit in probability of C £ {X,Y) as e — > + (of course, 
again under the condition of existence) . 

We recall that (cf. Protter [19]) a sequence of processes (H n ;n > 0) converges to a process H uniformly 
in probability if 

lim P ( sup \H n (t) - H(t)\ > a ) = for every a > 0. 
«->°o \te[o,T] J 

In [22] (Theorem 2.1) Russo and Vallois derived the ltd formula for the backward Russo-Vallois integral. 

Theorem 2.2. Let f 6 C 2 (M) and X be a continuous process admitting the generalized bracket, i.e., 
[X, X] exists in the sense of Definition 2.1. Then for every t e [0,T], the backward Russo-Vallois inte- 
gral J* f'(X(s))dX(s) exists and 

f(X(t)) = f(X(0)) + J" f'(X(s))dX( S ) \ J* F"(X(s))d[X, X](s), 

for every t > 0. 

We list some properties of Russo-Vallois integral, which will be used later in this paper. 

Proposition 2.3. (1). Lf X is a finite quadratic variation process (i.e., [X, X] exists and [X, X]t < +oo, 
P-a.s.) and Y is a zero quadratic variation process (i.e., [Y,Y] exists and equals to zero), then the mutual 
generalized bracket [X, Y] exists and vanishes, P-a.s. 

(2) . Lf X and Y have V-a.s. Holder continuous paths with order a and ft, respectively, such that a > 0, 
P > and a + (3 > 1, then [X, Y] = 0. 

(3) . We assume that X and Y are continuous and admit a mutual bracket. Then, for every continuous 
process {H(s) : s G [0,T]} 7 

f H(s)dC e (X,Y)(s) converges to [ H(s)d[X, Y](s). 
Jo Jo 

The following proposition, which can be found in Russo and Vallois [20] , states the relationship between 
the Young integral (see Young [23]) and the backward Russo-Vallois integral. 

Proposition 2.4. Let X,Y be two real processes with paths being P-a.s. in C a and C@ , respectively, with 
a > 0, j3 > and a + /3 > 1. Then the backward Russo-Vallois integral J Q YdX coincides with the Young 
integral J Q Yd^X. 
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3 A Generalized Ito Formula 



In this section we state a generalized Ito formula involving an Ito backward integral with respect to the 
Brownian motion W and the Russo-Vallois integral with respect to the fBm B. It will play an important 
role in our paper. It is noteworthy that this Ito formula corresponds to Lemma 1.3 in the paper of Pardoux 
and Peng [16] for the case of an fBm with Hurst parameter H = 1/2, i.e., when B is a Brownian motion. 

Theorem 3.1. Let a G C(H, [0, T]; R) be a process of the form 

a t = a + f Ads + [ jsl dW s , t G [0,T], 
Jo Jo 



where f3 and 7 are M-adapted processes and P{/ |/3 s |ds < 00} = 1 and P{/ |7s| 2 ds < +00} = 1, re 

i)y 



spectively. Suppose that F G C 2 (M x M). Then the Russo-Vallois integral L §p(a s , B s )dB s (defined as the 



1 ft 

■<-y oj 

surely, for all <t <T, 



uniform limit in probability of ^ J Q (B s — B s _ e )^(a s , B s )ds) exists forO <t<T, and it holds that, P-almost 



/"* OF f* dF f* dF 

F(a t ,B t )=F(a ,0) + ^ ^(a„ B s )(3 s ds + ^ — (a., B s ) 7s j dVF s + ^ — (a a ,B a )<LB a 



1 /"* c} 2 ^ 

-(a s ,B s )\j s \ 2 ds. 



2 7o 9a; 2 

Proof: Step 1. First we suppose Fe C 2 (K x M) (i.e., the function F is twice continuously diffcrentiable 
and has bounded derivatives of order less than or equal to two) and there is a positive constant C such that 
/ |/3 s |ds < C and J |7s| 2 ds < C. It is direct to check that 

1 r* 

F(a u B t )-F{a ,0) = lim- / (F(a s , B.) - F{a s _ £ , B s _ £ ))ds. 

e^O e Jo 

For simplicity we put a a , £ . s = a s - a(a s - a s _ £ ) and B a ,e, s = B s -a(B s - B s _ e ), for any a G [0, 1], s G [0,T], 
e > 0. We have 

F(a s ,B s )-F(a s _ £ ,B s _ £ ) 

=(a s - a s _ £ ) — (a s , B s ) + (B s - B s _ e ) — (a s ,B s ) - (a s - a s _ £ ) 2 J (a , e , s , £ ,e,«)(l - a)da 

f 1 a 2 F (3.2) 
- 2(a s - a s _ £ )(B s - B s _ e ) y ^^(a 0;£ , s , S ai£;S )(l - a)da 

-{B s -B s - £ f / — 1 -(a ,£, > ,fia,£, s )(l - a)da. 

jo cy 

By applying the stochastic Fubini theorem, we get that 
1 

-y (a a - a 5 _ £ ) — (a 5)J B s )ds 

= el (/ /3rdr + / > + dW ^) ^( a s,B s )ds 

I3 r — (a s ,B s )dsdr + - / / 7 r — (a s , B s )ds | dW r . 



Jr Cte e JO Jr ^ 

Since | J I f r+e ' )A * ^(a s ,B s )ds is % t -measurablc and converges to ^(a r ,B r ) when e — > 0, it follows that, 
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thanks to the continuity of 4^{a s , B s ) 



lim sup 

e ^°te[o,T] 



< lim sup 

£ ^°re[0,T] 



ft f I r(r+e)At Qp Qp \ 

J ^{ej r ^K,S s )d S -— (a r ,S r )Jdr 

/3 r |dr = 0, in probability. 
Thus, in virtue of the boundedness of by the Dominated Convergence Theorem, 

ft ( 1 f(r+e)M Qp Qp \ 



limE 

£^0 



I f(r+e)M Qp Qp 

— (a s ,B s )ds-—(a r ,B r ) 



sup 

te[o,T] 



< lim] 



sup 

r€[0,T] 



1 f(r+e)At Qp 



2 r* 




/ l7r| 2 dr 


J- T 


Jo 





0, P- a.s. 



(Recall that (± / r (r+e)A * §£ (a s , B s )ds - §§(a r , B r )) is H-adaptcd and W. - W T is an H-(backward) 

Brownian motion.) Thus, we get 



1 f OF f OF f OF 

lim - / (a s - a s - E ) — (a s ,B s )ds = / f3 r — (a r , B r )dr + / -f r — (a r , B r ) I dW r , (3.3) 
e^o e J dx J dx J dx 

uniformly in probability. We notice that the generalized bracket of a is the same as the classical one, i.e., 
[a, a] s = J Q S |7 r | 2 dr, s e [0, T]. We also have 



1 f f d F 

-J {a s -a s _ e ) 2 J -Q-^(a a ,e.s,B a ^^)(l - a)dads 

1 f* d 2 F 

~'2~e J ( a ^- a s-e) 2 - d —2(^s,B s )ds + A eit , te[0,T], 



(3.4) 



where 



A e , t = 



i r* 



1 (d 2 F 



d 2 F 



y dx2 -(a a ,e,s,B a!eiS ) - — {a Sl B s ) ) (1 - a)dads. 



Proposition 2.3 yields that ^ J ( a s — oe s - e ) 2 j^r(a s , B s )ds converges to J Q ^£(a s , B s )d[a, a] s and the con- 
tinuity of ^jt, a and B implies that A Sjt converges to zero. A similar argument shows that 



and the term 



If' f OF 

- / 2(a s - a s - e )(B s - B s - £ ) I (a ,e,a, B a , e ,s)(l - a)dads 

e Jo Jo oxoy 

1 I' Z" 1 d 2 F 

-J (B s - B s _ e ) 2 J (a a ,e,«,-B 0)e , s )(l - a)dads 



converge in probability, respectively, to 

32 



However, these both latter expressions are zero due to the fact that H e (1/2, 1) and Proposition 2.3 (Observe 
that the fBm has Holder continuous paths of any positive order less than H almost surely). 
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Combining the above results with (3.2), (3.3) and (3.4), we see that 

1 f t dF dF 

lim- J {B s -B s _ e ) — {a s , B s )ds =F(a t , B t ) - F(a o ,0) - J — {a s , B s )(3 s ds 

rt Qp ^ ft d 2 F 

-J ^K,S s ) 7s |dW s + - J -^{a s ,B s )\ ls \ 2 ds, 

uniformly in t G [0,T], in probability. Consequently, the integral J* ^(ot s , B s )dB s exists in Russo-Vallois' 
sense (Recall Definition 2.1) and we get the ltd formula (3.1) for F e C t 2 (1 x 1). 

Step 2. Now wc deal with the general case that P{/ T |/3 s |ds < oo} = P{Jq \"fs\ 2 ds < oo} = 1. For 
each n G N, wc define a sequence of H-stopping times by r n = sup{i < T : \/3 s \ds + J ( T |7 s | 2 ds > n} V 0, 
so we know that the processes {/?" := Pt\ Tn ,T] (*)> t G [0,T]} and {7" := j t l[ TnT ](t),t G [0,T]} satisfy 
Jq I ds < n and f Q T -7-™ | 2 ds < n, P-a.s. Furthermore, as n — > 00, r„ — >• 0, P-a.s. We consider the Ito 
formula for the process — a + J * /3™ds + J Q * 7™ I dW s ,t G [0, T]. Thanks to the result of Step 1, we have 
that, for every n G N, 

FK, B t ) =F(a , 0) + jT |^K, i? s )/3"d s + jf ^( a ^,B s )dB s + J* ^-( a n s ,B s ) ls | dW s 



-I 

2 Jo 



1 d 2 F 



9x 2 



(a?,B«)|7a|d«- 



Since a n converges to a uniformly in probability on [0,T], by letting n — ► 00 in the above equation, we 
deduce that 

lim / — «,B s )dB 5 = lim lim / - — «, B S )(B S - B s _ £ )d,s = lim / - — (a 8 , B S ){B S - B s _ £ )ds 
n^ooJ Q dy n^oos^0j e dy e^oJ Q e dy 

exists, it is the Russo-Vallois integral J Q §^(a s , B s )dB s and it equals to 

F(at,B t )-F(ao,0)- J -^(a s ,B s )f3 s ds - J -(a 8 ,B s ) 7 4d^ + -| ^ (a s , B s )| 7s | 2 d S . 

Step 3. Finally we consider the case F G C 2 (R x R). We let {(^jvIatgn be a sequence of infinitely 
diffcrcntiable functions with compact support such that ip^ix) = x for {(xi,X2) '■ max(|xi|, X2I) < N},N G 
N. We set F N (x) = F(<p N (x)), so that F N (x) G C b 2 (R x R) for every N > 0. We notice that F N (a.,B.) and 
F(a.,B.) coincide on the set SIn = G f2 : sup 0<s<t a s < iV,sup 0<s<t £? s < V} and that = Uat>i ^n- 
Due to Step 1 and Step 2, for every N, we have 

ft €\t^ ft c\ Tp rt c\r? 

F N (a t ,B t )=F N (a ,0) + J -^(a s , B s )f3 s ds + jf -^(a 8) B s )dB s + J -JL(a s ,B s ) ls idW s 
-\f o ^{a s ,B s )\ ls \ 2 ds, te[0,T\. 

Therefore, for every N, it holds on that 

f* dF f* dF f* dF 

F(a t ,B t )=F(a ,0) + j —{a s , B s )/3 s ds + j ^(a„ B s )dB s + j —(a s , B s ) ls I dW s 

Finally, by letting AT tend to +00, we get the wished result. The proof is complete. ■ 
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4 Doss-Sussman Transformation of Fractional Backward Doubly 
Stochastic Differential Equations 

In what follows, we use the following hypotheses: 

(HI) The functions a : W l -> R nxd and b : R n -> R™ arc Lipschitz continuous. 

(H2) The function / : [0, T] x M™ x M x R d -> E is Lipschitz in (x, y, z) G K" x M x R d -> R with 
|/(t, 0, 0, 0)| < C uniformly in t G [0,T], the function 3 : R™ -> K belongs to C fc 3 (M) and the function $ is 
bounded. 

We hx an arbitrary t G [0,T] G K+. Let (Xl' x ) < s < t be the unique solution of the following stochastic 
differential equation: 

f dX s *>* = -&(* a t>x )d* - tr(X^) | dW s , s G [0, i] , 

\ = x. (4 ' 1} 

Here the stochastic integral J Q ■ 1 dW s is again understood as the backward ltd one. The condition (HI) 
guarantees the existence and uniqueness of the solution (^' x )o<s<t in M 2 (¥ w , [0,T];M"). Our aim is to 
study the following backward doubly stochastic differential equation: 

Jo Jo Jo 

We emphasiz that the integral with respect to the fBm B is interpreted in the Russo-Vallois sense, while 
the integral with respect to the Brownian motion W is the ltd backward one. If B is a standard Brownian 
motion, equation (4.2) coincides with the BDSDE which was first studied by Pardoux and Peng [16] in 1994 
(apart of a time inversion). 

Before we investigate the BDSDE (4.2), we first give the definition of its solution. 

Definition 4.1. A solution of equation (4.2) is a couple of processes {Ul' x , V s *' x ) se [o l t] such that: 

1) . {U^X' x )sem G ^r(l) x H 2 t {R d ); 

2) . The Russo-Vallois integral J Q g(U* ,x )dB r is well defined on [0,t]; 

3) . Equation (4.2) holds W-a.s. 

Unlike the classical case, the lack of the semimartingale property of the fBm B gives an extra difficulty 
in solving BDSDE (4.2) directly. However, the work of Buckdahn and Ma [5] indicates another possibility to 
investigate this equation: by using the Doss-Sussman transformation. Let us develop the idea. We denote 
by n the stochastic flow which is the unique solution of the following stochastic differential equation: 



v(t,y) = y+ [ g(v(s,y))dB s , t e [0,T], (4.3) 
Jo 

where the integral is interpreted in the sense of Russo-Vallois. The solution of such a stochastic differential 
equation can be written as r](t,y) = a(y,Bt) via Doss transformation (see, for example, Zahlc [25]), where 
cx(y, z) is the solution of the ordinary differential equation 

^(y,z)=g(a(y,z)), z G K, (4 4) 



a 



M) 



By the classical PDE theory we know that, for every z G R, the mapping y h-> a(y, z) is a diffeomorphism 
over R and (y, z) i-> a(y, z) is C 2 . In particular, we can define the y-inverse of a(y, z) and we denote it by 
h(y, z), such that we have a(h(y, z)) — y, (y, z) G R x M. Hence, it follows that 

^(h(y,z),z)^(y,z) = 1 and ^(h(y,z),z) + ^(h(y,z),z)^(y,z) = 0. 
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Therefore, 



As a direct consequence we have that also T)(t,-) — a(-,B t ) : R — > R is a difFeomorphism and, thus, we can 
define S(t,y) := r)(t, = Hy> B t),(t,y) G [0,T] x R. Moreover, by the Ito formula (Theorem 3.1), we 

have 

d£(t,i/) = d%,B t ) = — h{y,B t )dB t = -—£{t,y)g(y)&B u t e [0,T], 
i.e., the process £ satisfies the following equation: 



jf ^£(s,y)g(y)dB s , t e [0,T]. 



(4.5) 



Furthermore, we have the following estimates for 77 and 5. 



Lemma 4.2. There exists a constant C > depending only on the bound of g and its partial derivatives 
such that for £ = 77, £, it ZioZds £/iai, P-a.s., for all (t, y), 

m,y)\<\y\+C\B t \, exp{-C|S t |}< M <exp{C\B t \}, 



<),i- 



<exp{C|i? t |}, 



<exp{C|i? t |}. 



Proof: The first three estimates are similar to those in Buckdahn and Ma [5]. So we only prove the last one. 
For this end we define 7(6*, y, z) = a(y, 6z), for (9, y, z) e [0, 1] x R x R™. It follows from (4.4) that 

f 6 

l{B,y,z) = y+ I (g(-y(r,y,z)),z)dr. 
Jo 

By differentiating this latter equation, we obtain, 



£fa{6,y,z)= fymv,z))z($(6,y,zj) +3$fc(>y(6,y,z))z%(6,y,z)${e,y,z) 

0(0,^) = o, 



(4.6) 



and from the variation of parameter formula it follows that 



~J^' :L !I- : ' = J q CXp {/ ^' : " ' :, '--' / - ::::(ir )f "^ ( ~ :::: ! :i 



dy 



^9, 



d 2 l, 



Thus, by using the first three estimates of this Lemma, we get 

d 3 



Q3 



<exp{C|B t |}. 



Hence we have completed the proof. 

Lemma 4.2 plays an important role in the rest of the paper thanks to the following lemma. 



Lemma 4.3. For any C G R, we have 



E 



exp < C sup \B S 

I «€[0,T] 



< OO. 
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Proof: The proof is similar as Lemma 2.4 in [11] (and even easier), so we omit it. 

We denote by fl' the subspace of 0' such that 0' = jo/ G 0' : sup sG [ T ] < oo 
P'(Q') = 1. 

We let (Y*' x , Z*- x ) be the unique solution of the following BSDE: 



- <&(^' x ) + / /(r, X*'*, y r *'*, Z*>*)d 



where 



f(t,x,y,z) 



I | t,x,rj(t,y), ^-rj(t,y)z ) -U: 



idW r 



d 2 



is clear that 



(4.7) 



This BSDE, studied over fi = fi' x Q" and driven by the Brownian motion W r (uj) — W r (w") = w"(r), 
r G [0,T], can be interpreted as an w'-pathwise BSDE, i.e., as a BSDE over f2", considered for every fixed 
lj' G 0'. However, subtleties of measurability make us preferring to consider the BSDE over Q, with respect 
to the filtration H. We point out that the coefficient / has a quadratic growth in z, while the terminal value 
is bounded. BSDEs of this type have been studied by Kobylanski [12]. We state an existence and uniqueness 
result for this kind of BSDE, but with a slight adaptation to our framework. For this we consider a driving 
coefficient G satisfying the following assumptions: 

(H3) The coefficient G : ft x [0, T] x R x R d is measurable, for every fixed (y, z), progressively measurable, 
with respect to the backward filtration H and G is continuous in (t, y, z); 



(H4) There exists some real-valued -measurable random variable K : ft' 
K(l + \z\ 2 ). 



such that \G(t,y, z)\ < 



(H5) There exist real-valued -measurable random variables C > 0, e > 0, and ® £>([0, Im- 



measurable functions k, l £ : fi' x [0, T] — >• M such that 

<k(t) + C\z\, for all {t,y,z), P-a.s., 



dG, 



dG 

— (t,y,z)<l e (t) + e\z\ 2 , for all (t,y,z), P-a.s. 

Remark 4.4. Due to the Lemmata 4-2 and 4-3, the function f in the equation (4.7) satisfies (H3) — (H5). 
In particular, K = exp{Csup ig [ T ] \B t \} for C G M + appropriately chosen. 

Adapting the results by Kobylanski [12] (Theorem 2.3 and Theorem 2.6), we can state the following: 

Theorem 4.5. Let G be a driver such that hypotheses (H3)-(H5) hold and let £ be a real-valued 'HQ- 
measurable random variable, which is bounded by a real-valued -measurable random variable. Then there 
exists a unique solution (Y, Z) G Uf{M) x H^(R d ) of BSDE 



Y t = £+ f G{s,Y s ,Z s )ds- f Z s idW s , t£[0,T]. 
Jo Jo 



(4.8) 



Moreover, there exists a real-valued -measurable random variable C depending only on esssup[ T ] xr2 /, 
\Y t (uj', to")\ andK, such that 

fT 



E 



\Z s \ 2 ds\rf 



a.s. 
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Remark 4.6. The conditional expectation E [-IJ 7 ^] is here understood in the generalized sense: if £ is a 
nonnegative TLo -measurable random variable, 

E[£|7f] := lim t E[£ An\F^](< oo) 

n— t-oo 

is a well defined -measurable random variable. If £ is noi nonnegative we decompose £ = £ + — = 
max{£,0},£- = -min{£,0} and we put E[£| Jf] = E[£+|7^] -E[£-|J"^] on {min{E[£+|.F^],E[£-|7f ]} < 
oo}. 

Proof of Theorem 4.5: We observe that the Brownian motion W possesses the (backward) martingale 
representation with respect to the backward filtration H, i.e., given an 'Ho-measurable random variable £ 
such that E [^ 2 |^ r y ] < oo, P-a.s., there exists a unique process 7 £ V.^(M. d ) such that 

£ = E[f I ] - / jri dW r , P - a.s. 
Jo 

This martingale representation property allows to show the existence and uniqueness of a solution of (4.8) 
when G is of linear growth and Lipschitz in (y,z). Combining this with the approach by Kobylanski [12] 
allows to obtain the result stated in Theorem 4.5. ■ 



By using Theorem 4.5, we are now able to characterize more precisely the solution of BSDE (4.7). 

Theorem 4.7. Under our standard assumptions on the coefficients a, b, / and $, BSDE (4.7) admits a 
unique solution (Y t ' x 7 Z t ' x ) in 7i%°(M) x Tl 2 (R d ). Moreover, there exists a positive increasing process 9 £ 
L°(H,R) such that \Y*' X \ < 9 S , E [J Q T \Z^ x \ 2 As\H T ] < exp{exp{Csup se[M \B S \}}, F-a.s., for all U-stopping 
times t (0 < r < t), where C is a constant chosen in an adequate way. Furthermore, the process (Y t,x , Z t,x ) 
is G -adapted. 

Proof: Due to Theorem 4.5, equation (4.7) has a unique solution Z Lx ) in Hf°{R) x W. 2 (R d ). 

Step 1. In order to give the estimates of {Y t,x , Z l - X ), we proceed as in Lemma 5.3 in [5]. In particular, 
we can show that there exists an increasing positive process 9 £ L°(¥ B , [0, T]) such that P-a.s., \Y*' X \ < 
9 S ,0 < s < t < T. This process (9 S ) can be chosen as the solution of the following ordinary differential 
equation 

r\9 

-±=exp{C sup \B r \}(l + 6 s ); 6(0) = 

OS 0<r<t 

for some suitably chosen real constant C, i.e., 

9 S = (|$(xS ,x )| + l)exp|exp|C o SupjB r ||*| - 1. 

Indeed, for M > 0, let (pmiy) be a C°° function such that < ipu < 1, Vuiy) = 1 for \y\ < M and 
1 -Pm(v) = for 1 2/ 1 > M + l. Defining a new function f M by f M (t, x, y, z) = f(t, x, y, z)<Pm(u) we see that the 
function f M also satisfies conditions (H3)-(H5). According to Theorem 4.5, there exists a unique solution 
(Y M ' t < x ,Z M ' t - x ) of equation (4.7) with / being replaced with f M . The stability result in Kobylanski [12] 
shows that, when M — > +00, there exists a subsequence of Y M ^' X converging to Y l - X uniformly in probability. 
Therefore, following Buckdahn and Ma's approach and slightly adapted, we only need to prove that Y^ r ' t ' x 
is uniformly bounded by 9 S . We apply the Tanaka formula to \Y M ' t ' x \ to get that 

|yM, t ,x| = f sgn (Y r M ^ x ) f M (r,Xl' x ,Y r M ^ x ,Z^^ x )Ar- f sgn (Y r M '*' x ) Zf^ x 1 dVK r +L s -L , 

Jo Jo 

for a local-time-like process 1/ such that L t — and L s = J g l^ y M,t,x =0 jdL r . Now we define a new function 
ijj(y) by letting ?/>(?/) = e 2Ky — 1 — 2ifj/ — 2K 2 y 2 , for y > 0, and = for y < 0, where K is the bound 
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in Remark 4.4. Then we apply the Ito formula to ^>(|Y S M '*' X | — S ) and get 
^(\Y S M ^\-6 S ) 

= f i,' (\Y r M ^ x \ -e r )sgn(Y r M ^) f/ M (r,X^,y r M ^,Z r M ^)-exp{C sup |B U |}(1 + r )) dr 
Jo V o<«<t / 

- f 4>' (\Y r M ^ x \ - 6 r ) sgn (Y r M ^ x ) Z r M ^ x | dW r + f {\Y r M ^ x \ - 9 r ) dL r 
Jo Jo 

-\J^" (\Y r M ^ x \ - 6 r ) \Z^ x \ 2 dr 
The property of ip shows that f° V' (|F r M '*' x | - 6 r ) dL r = f° tp' (-6 r ) dL r = 0. We also have 

/V (|n M ''i - Or) sgn(Y r M '^) (f M (r,Xl' x ,Y r M ^,Z^ x ) - exp{C sup \B u \}(l + 8 r )) dr 

JO \ 0<u<t J 

< [' V' (\Y r M ^ x \ - r ) (K {^{Y r M ^ x ) \Y r M ^ x \ - e r ) + K\Z^ X \ 2 ) dr. 
Jo 

Since ip" — 2Kip' > 0, we get from equation (4.9) that 

V> (|y s M '*' x | - S ) <K / S V (Ir^l - T ) (<p(Y r M '*' x ) \Y r M ^ x \ - r ) dr 
Jo 



(4.9) 



Jo 



M,t,x\ 



B r ) sgn (Y r ' t,x ) Z^ x | dW, 



Thus, we deduce that 

E[V(|n M '*' x | -0 s )\n a ] <e \ K [%/ (\Y r M ^ x \ -6 r ) (ip(Y r M ^ x )\Y r M ^ x \ -9 r )dr 

L Jo 

From the definition of tp we get that "4>'(y) — 2K(tp(y) + 2K 2 y 2 ). Hence, we have 
E[i>(\Y s M ^ x \-6 s )\H s ] 

<E [ [ S 2K 2 i>(\Y r M > t ' x \ -6 r ) (y{Y r M ^ x )\Y r M ^ x \ ~ 9 r ) +4K 4 (\Y r M ^ x \ - 8 r ) 3 dr 
. Jo 

There also exists a K such that y 3 < Kip(y), for all j/gl. Consequently, we get that 

E[i,(\Y s M ^ x \-e s )\n s ] <2K 2 {M+\\6\\ 00 , m +2K 2 K) [' E [V (|F r M '*' x | - 9 r ) \ H s ] dr. 

Jo 



Finally, the Gronwall inequality shows that ip (\YJ M ' t ' x \ — S ) = 0, for any s G [0,i], P-a.s. Therefore, 
|yM,t,x| < ^ for any g £ [ ,t], P-a.s. 

Step 2. We apply the Ito formula to e ay =' x , with a being a real-valued J 7 ^ -measurable random variable 
to be determined later, and we obtain 

e «e =e «*(4') + [' ae aY 'f(r,Xt> x , Y?>*, f \a 2 e aY r" \Z^ x \ 2 dr - f ae aY r x Z^ x | dW r 

Jo Jo 2 J 



< e »*(X*-) + 



1 



^^ + |a|X) e ay -' x |y r t ' :E | 2 dr+ / lal^e^'dr - / ae ar '" Zf.' x | dW r . 



Jo 



o aY*' X ryt.X 



Let C(w') := esssup[Q t ] x q» |l^* ,:c (a/, u/') | (< +oo,a/ G 0'). Hence, by taking the conditional expectations 
E[-|%t] on both sides, we deduce that for any H-stopping times r £ [0,T], 



\Zp x \ 2 dr 



Qa 2 - \a\Kj e-^E J 
Jo 



< [-a 2 -\a\K )E 



|Z r | 2 dr 



7^ T 
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We can choose a 



\a 2 — \a\K — AK 2 and we get, keeping in mind that here K is a random 
variable bounded by exp{C sup se [ t ] \B S \} (C G R + is a real constant, see Remark 4.4), 



E 


f r \zr\ 2 dr 






.Jo 





< {2 + AKH \ ^ < exp I exp { C sup \B S 



4K 2 



s€[0,t] 



(4.10) 



Step 3. Let us show that the process (Y f ' x , Z f ' x ) is not only H- but also G-adapted. For this we consider 
for an arbitrarily given r G [0,t] equation (4.7) over the time interval [0,r]: 

Y*' x = $(X*' X ) + f f (r, X' r <*, Y r *' x , Z?) dr - f Z\> x | dW r , s G [0, r] . (4.11) 
Jo Jo 

Let ri T t := T^t v J~? i * e [0i T ]- Then H T = {'H^}te[o.T] is a backward Brownian filtration enlarged by a a- 
algebra generated by the fBm B, which is independent of the Brownian filtration. Thus, with respect to HT, 
the Brownian motion W has the martingale representation property. Since f(r, x, y, z) is Qt- and, hence, also 
"H^-measurable, dr a.e. on [0, r], it follows from the classical BSDE theory (or Theorem 4.5) that BSDE (4.11) 
admits a unique solution (Y*^, Z*^) G ft~(R) x H 2 (R d )- On the other hand, also (Y*> X > T , Z^^)^ _ r] 
is a solution of (4.7). Hence, (Y*' X ,Z*' X ) = (Y*' X ' T , Z*' x ' r ), dr a.e., for i < r. Consequently, (F r *> x , Z*- x ) 'is 
"H^- measurable, dr a.e., for r < r. Therefore, letting t jt we can deduce from the right continuity of the 
filtration F B that (Y*< x , Z^ x ) is G-adapted. ■ 

Remark 4.8. We remind the reader that the bound we get in (4.10) is only P-o.s. finite, but not square- 
integrable. As a matter of fact, it is hard to prove directly that Z l - X is a square-integrable process, which 
constitutes the main reason that we use instead the space 'H 2 (M. d ). Hence, the major difference between our 
work and Buckdahn and Ma [5] is: In the classical case, a priori we can solve the BDSDE in the first step 
to get the square integrability of Z , but in the fractional case, there is not a direct way to solve the BDSDE. 

Now we are ready to give the main result of this section by linking the BSDEs (4.2) and (4.7) with the 
help of the Doss-Sussman transformation. 

Theorem 4.9. Let us define a new pair of processes (£7*' x , V f ' x ) by 



d 



where (T*' 1 ,^ 1 ; 
BDSDE (4.2). 



U t s > x = V (s,Y s t > x ), 
is the solution of BSDE (4.7). Then (U^ x , V f ' x ) G U 



xH 2 



is the solution of 



Remark 4.10. The above theorem can be considered as a counterpart of Theorem 3.9 in Jing and Leon [11] 
for the semi-linear case when H < 1/2. However, since we use here a different Hurst parameter H and a 
different type of stochastic integral with respect to fBm B, the above theorem obviously does not cover the 
result in [11]. 

Proof of Theorem 4.9: The fact that (t/*' x , V t,x ) G %f(R) x H 2 (R d ) follows directly from Theorem 4.7 and 
Lemma 4.2. In order to prove the remaining part of the theorem, we just have to apply the ltd formula 
(Theorem 3.1) to a(Y*' x ,B s ), noticing a(Yf> x ,B 3 ) = r](s,Y*' x ), to obtain that, for (s,x) G [0,t] x R d , the 
Russso-Vallois integral J" s g(Up x )dB r exists and 



1 

+ 2 tr 



d_ 

>y 

d 2 



—n(r Y l ' x ) 



(Z^f^virXnZ^ 



dr - Z^ x I dW r 



tr 



(Z^f^rjirX^Zy- 



+ 



f 

Jo 



d_ 



a{Y*' x ,B r )dB r . 



dr 



(4.12) 
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Consequently, 

Ut' x = + / f(r,X t r ' x ,U t r ' x ,V*> x )dr + f g(U}' x )dB r - f V^ x idW r . (4.13) 

Jo Jo Jo 

The proof is complete. ■ 
Now we close this section by a simple example to illustrate the idea of the above procedure. 

Example 4.11. Let us consider the linear case ( for simplicity of notations we omit the superscript (t,x)): 

f dU s = U s dB s + f(U s ,V s )ds + V S I dW s , s G [0,i\; 

{ u = (4J4) 

It is elementary to show that a{y,z) = ye z . Hence, the solution (U, V) of equation (4.14) is given by 
U s = Y s e Bs and V s — Z s e Bs , where (Y, Z) is the solution of 

f dY s = f(Y s , Z s )ds + Z S { dW s , s G [0, t}; 

I Y = (415) 

and the function f is defined by f(y, x) = e~ Bt f(ye Bt , ze Bt ). Obviously f is Lipschitz in (y, z) as far as f 
is Lipschitz. By following a classical Malliavin calculus method {see, e.g., Pardoux and Peng [16]), we can 
get that Z is P-a.s. uniformly bounded and, thus, the process {Y s ,s <G [0,t]} is (1/2 — e)-Hdlder continuous 
in s, for all e G (0, 1/2). Therefore, we have, for every r, s G [0,t], 

\U r -U s \ = \Y r e B - -Y s e B °\ < ||F|U|e Sr - e B ° \ + el B =l \Y r - Y s \ 

<cxp{ sup |B u |l||y|| 0O C(u/)|r--s| a +exp{ sup \B U \ \ C{u')\r - s| 1/2 " e , for a.a.w' G Q! , 

where we can choose a to be in (1/2, H). That is to say, we can choose < e < a — 1/2 and get that the 
process {U s ,s G [0,t]} has a-Holder continuous paths. Hence, instead of using the Russo-Vallois integral 
with respect to the fractional Brownian motion in equation (4.14) we can use the classical Young integral. 
Furthermore, thanks to Proposition 2.3, these two integrals coincide. 

Remark 4.12. In the latter example, the Young integral and the Russo-Vallois integral coincide. Un- 
fortunately, the Holder continuity seems to be very hard to deduce in the general, nonlinear case. As a 
consequence, we have to work with the more general Russo-Vallois integral. 



5 Associated Stochastic Partial Differential Equations 

In this section we discuss briefly the relationship between an associated SPDE and a PDE with stochastic 
coefhents. For simplicity, we only show the relationship for the case of classical solutions. For a complete 
discussion of the case of (stochastic) viscosity solutions, one can proceed by adapting the approaches in 
Buckdahn and Ma [5] and Jing and Leon [11]. 

Let ££ be the second order elliptic differential operator: 

1 ™ rfi n f) 

1,3 = 1 i= 1 

which means that it is the infinitesimal generator of the Markovian process {X*' x ,s G [0,t]} defined by 
equation (4.1). 

Our aim is to study the following semilinear SPDE driven by the fBm B: 

f dv{t,x) = [Sfv(t,x) - f(t,x,v(t,x),<j(x) T J^v(t,x))] dt + g(v(t,x))dB t , (t, x) G (0, T) x R"; , . 
\ v(0,x) = $(x), x G R n . y ' 
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In the case that B is a Brownian motion (an fBm with Hurst parameter H — 1/2), it is well known (see, 
Pardoux and Peng [16], Buckdahn and Ma [5]) that the random field v(t, x) :— U\ ,x solves (5.1) (in the 
viscosity sense if the coefficients are Lipschitz, and in the classical sense if the coefficients are C|), where 
U t,x is the solution of BDSDE (4.2). Thus, it is natural to raise the following question: Can we also solve 
the SPDE (5.1) driven by the fBm B, by studying the properties of the solution of its associated BDSDE 
(4.2)? The answer is positive. Indeed, by applying the Doss-Sussman transformation, we will show that the 
PDE: 



du(t,x) = (&u(t,x) - f{t,x,u{t 1 x) 1 cj{x) T -^u{t,x))^ dt, {t,x)€(0,T) 



u(0,x) = ${x), 
is transformed into SPDE (5.1), where we recall that 



x 6 



(5.2) 



d 

f(t,x,y,z) = [ Q^v(t,y) 



d 1 

fit, x, r](t, y), g^v{t, y)z) + -TV 



d 2 



We also observe that following a similar argument as in Kobylanski [12], under some smoothness assumptions, 
u(t,x) := F t M , (t,x) G [0,T] x K" is the solution of equation (5.2), where Y f ' x is the solution of BSDE (4.7). 
First we give the definition for the classical solutions of equations (5.1) and (5.2). 

Definition 5.1. We say a stochastic field w : Q' x [0, T] x M™ — > M is a classical solution of equation (5.1) 
(resp., (5.2)), if w G 2 and satisfies equation (5.1) (resp., (5.2)). 

We have the following proposition. 

Proposition 5.2. Suppose that u is a classical solution of equation (5.2), then u{t,x) = n(t,u(t,x)) = 
a(u(t,x), Bf) is a classical solution of SPDE (5.1). The converse holds also true: Every classical solution it 
of equation (5.1) defines a classical solution u(t,x) = £(t,u(t,x)) of equation (5.2). 

Proof: The claim that u(t, x) £ C ¥ ' B follows from the regularity property of the functions a and u and the 
fact that u is F B -adapted. Moreover, we first observe that u(0, x) = rj(Q,u(0,x)) = u(0,x) = $(x), and we 
apply the ltd formula to tt(t, x) to obtain: 



au(t,x) =-^-a 
dy 



d 



(u(t,x),B t ) (^u{t,x) - / ^t,x,u(t,x),o-(x) T ^-u(t,x)^ 



dy 



c) 



d 

dt+ — a(u(t,x),B t )dB t 



d 



—a(u(t, x),B t )Jfu(t, x) — f I t, x, r)(t, u(t, x)), -^-f}{t, u(t, x))a(x) —u(t, x) 



-Tr 



dy 



d d 2 d 

ct(x) t — u(t, x)—ri{t, u(t, x))—u(t, x) T a(x) 



dx 



dt+ — a{u{t,x),B t )dB t . 
az 



(5.3) 



We notice that 







d 



—u(t,x) = —[a(u(t,x),B t )] = —a (u(t,x), B t )—u(t,x), 



dx 



dx 



d 



dy 



d 



dx 



d 2 d ( d 2 \ ( d \ T ( d \ d 2 

dx 2 ^ 7 x) = x) \dy 2 ' a ) ^ Bt) \d^ u ^ 7 x) ) + \dy~ a ) ^ x) ' ^dx 2 ^' x) - 

Thus, by substituting the above relations into (5.3), we can simplify the equation (5.3) and we get 



du(t, x) 



'/ u ( / . .i ) / | t, x, u(t, x), a(x) T -^-ii(t, x) 



dt + g(u(x,t))dB t . 



(5.4) 



2 CfB i s the space of F s -adapted processes u(t, x) which arc continuous in t and C 2 in x. 
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Consequently, u is a solution of SPDE (5.1). The proof of the converse direction is analogous. The proof is 
complete. ■ 
In order to summarize, we have constructed a solution of SPDE (5.1) with the help of the fractional 
BDSDE (4.2), by passing through the quadratic BSDE (4.7) and the associated PDE (5.2) with random 
coefficients: 

fract. BDSDE (4.2) T/l ^ 49 quadratic BSDE (4.7) 

fract. SPDE (5.1) Pr ^- 2 PDE ( 5 . 2 ) (with random coefficients). 

Proposition 5.2 shows that, in the case of a classical solution, the Doss-Sussman transformation establishes 
a link between PDE (5.2) and SPDE (5.1). This motivates us to give the following definition of the stochastic 
viscosity solution. For the cases H e (0, 1/2], the reader is referred to Buckdahn and Ma [5], Jing and Leon 
[11]. For the classical definition of the viscosity solution, we refer to Crandall ct al. [6]. 

Definition 5.3. A continuous random field u : fi' x [0, T] x E™ — > E is called a (stochastic) viscosity solution 
of equation (5.1) if and only if u(t,x) = £(t,u(t,x)), (t,x) 6 [0,T] x M™ is the viscosity solution of equation 
(5.2). 

By following a similar argument as that developed in the proof of Theorem 4.9 in Jing and Leon [11], 
our preceding discussion leads to the following theorem: 

Theorem 5.4. The stochastic field u : 0' x [0, T] x E™ — > E defined by u(t, x) = a(Y^ x , B t ) is a (stochastic) 
viscosity solution of SPDE (5.1). 
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